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Abstract 

Heterogeneous networks are a flexible deployment model that rely on low power nodes to improve the 
user broadband experience in a cost effective manner. Femtocells are an integral part of heterogeneous networks, 
whose main purpose is to improve the indoor capacity. When restricting access to home users, femtocells cause a 
substantial interference problem that cannot be mitigated through coordination with the macrocell base station. In 
this paper, we analyze multiple antenna communication on the downlink of a macrocell network, with femtocell 
overlay. We evaluate the feasibility of limited feedback beamforming given delay on the feedback channel, 
quantization error and uncoordinated interference from the femtocells. We model the femtocell spatial distribution 
as a Poisson point process and the temporal correlation of the channel according to a Gauss-Markov model. We 
derive the probability of outage at the macrocell users as a function of the temporal correlation, the femtocell 
density, and the feedback rate. We propose rate backoff to maximize the average achievable rate in the network. 
Simulation results show that limited feedback beamforming is a viable solution for femtocell networks despite the 
CSI inaccuracy and the interference. They illustrate how properly designed rate backoff improves the achievable 
rate of the macrocell system. 
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I. Introduction 

The best way to increase the capacity of cellular systems is by decreasing the cell size [[TJ-[[3j. 
Unfortunately, relying on cell splitting in traditional network deployment approaches is complex in high 
density environments. Heterogeneous networks provide a more flexible deployment model that relies 
on base stations of diverse properties to improve the user broadband experience in a ubiquitous and 
cost effective manner |4j. They consist of macro base stations overlaid with low power nodes such 
as pico base stations, distributed antennas, relays, and (of primary interest in this paper) femtocells. 
Unlike distributed antennas and relays, femtocells are user deployed and unplanned by the network 
operators. They are short range, operate on licensed spectrum and are usually restricted to serve home 
or enterprise users only. Unfortunately, restricting access to the femtocells to subscribed home users 
causes significant interference between adjacent femtocells, as well as cross tier interference between 
the femtocells and the macrocells [[3j, (5J. Moreover, interference management in femtocell networks 
does not conform with traditional cellular networks [[6j, |7J, as it cannot rely on coordination between 
the macrocells and the femtocells and needs to be decentralized. 

As interest in multiple input multiple output (MIMO) communication has grown, upcoming cellular 
standards have embraced multiple antennas on the base stations and the mobile users to increase data 
rates and improve performance of the radio link (5J. Multiple antennas are being considered at the 
femtocells for coverage improvement and interference management [[3j, [|9j. One of the most practical 
transmission strategies for realizing the benefits of MIMO systems is limited feedback linear precoding 



1 10 1 . Limited feedback suffers, however, from a CSI inaccuracy due to feedback delay and quantization 



error [11|, [12|. Its performance further degrades in cellular systems, when the interference is not 
coordinated, and the transmitter does not have instantaneous knowledge about the interference at the 
receiver JT3J. In this case, the transmitted rate does not take into account the accurate interference 
information, resulting in an increased probability of outage. In this paper, we consider limited feedback 
beamforming on the downlink of a heterogeneous network with private access femtocells, and we 
evaluate its performance, given the cross tier interference in the network, and the feedback delay. 

A. Prior Work 

Prior research on femtocell networks has focused on decentralized and self-optimization strategies for 
interference management as well as femtocell access control [14|-[19|. The feasibility of deploying 



femtocells in the same frequency bands as the existing macrocell, and methods for femtocell power 
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auto-configuration and public access were investigated in |2j, p"5| . In flT6| , the authors analyzed the 
effect of the additional handovers caused by femtocells on the dropped call rate of the macrocell 
users. Interference avoidance techniques such as self-optimization and dynamic frequency planning 
were suggested in fT8| to improve the capacity and decrease interference for an OFDMA femtocell. 
The performance of the femtocell access control strategy was investigated in JI9), with multi-cell and 



single-cell processing. Femtocell access strategy was discussed in [17|, where it was shown that for 
OFDMA femtocells, the femtocell access control should be adapted to the cellular user density. 

Thus far there has been limited work on the application of multiple antennas to heterogeneous 
networks with femtocells. The benefits of using multiple antennas at the femtocells and the macro 
base stations, in terms of coverage and spatial reuse, were evaluated in [9j. The analysis, based on 
zero forcing precoding, assumed full CSI at the transmitter, and suggested that single user transmission 
outperforms multi user transmission on the downlink of heterogeneous networks. Precoding techniques 
at the femtocells, based on having perfect or quantized CSI at the transmitter, were investigated in 



1 20 1, [21 1 to reduce the interference, perform local coordination, or analyze the best transmission 
strategy for the femtocells. The work considered two adjacent femtocell base stations only, and failed to 
capture the randomness in the deployment of the femtocells, as well as the feedback delay. Furthermore, 
evaluating the effect of cross tier interference on the performance of the cellular network borrows from 



the literature on transmission capacity for ad hoc networks with multiple antennas [22|-[26|, where 
the outage probability and transmission capacity for multiple antenna diversity and spatial multiplexing 
techniques were derived. An expression for the asymptotic spectral efficiency in an interference limited 
regime was derived in pTJ. The analysis considered interference from spatially distributed nodes with 
CSI at the transmitter, and concluded that CSI increases the spectral efficiency. 

To study the feasibility of limited feedback beamforming on the downlink of a heterogeneous network, 
the analysis should take into account the effect of feedback delay and quantization error. Feedback delay 
arises from signal processing algorithms, propagation, and channel access protocols. The achievable rate 
and bit error rate performances of MIMO systems with feedback delay have been investigated in pT|, 



[12], [|28|-[[3T|. The delay has been found to reduce the achievable throughput [ |29| , pi] , and to cause 
interference between spatial data streams [12]. Delay has been shown in [[30| to incur a significant 



performance loss on multi user MIMO systems, due to residual inter user interference, and in [32| to 
limit the performance of joint precoding and scheduling for the MIMO broadcast channel. The joint 



effect of delay and other cell interference assuming one strong interferer has been shown in [13| to 
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increase the decay in the achievable rate using a Markov chain model. In this paper, to investigate the 
applicability of limited feedback methods, we consider the joint effect of feedback delay and cross tier 
interference. This has not been investigated in the literature for randomly distributed interferers. 



B. Contributions 

In this paper, we consider limited feedback beamforming on the downlink of a heterogeneous network, 
with a femtocell overlay. We compute the probability of successful reception and the achievable rate 
for a macrocell user in the presence of cross tier interference from the femtocells. We account for the 
joint effect of feedback delay, quantization error and cross tier interference as follows. The propagation 
information that reaches the base station, being quantized and delayed, causes a CSI mismatch between 
the transmitter and the receiver. The mismatch translates into packet outage when the transmitted rate 



exceeds the instantaneous mutual information of the channel [13], [33 1, [34]. When the interference 



is not coordinated, the transmitter is oblivious to the instantaneous interference at the receiver, and 
modulates its information at a rate corresponding to the estimated signal to interference ratio (SIR) 
rather than the actual SIR, hence increasing the probability of outage. 

We consider a narrowband channel model with a Gauss Markov temporal correlation and a distance 
dependent path loss model. Since we compute the average rates as a function of the distance and we 
average over all locations inside the macrocell, we neglect large scale fading. The base stations have 
multiple antennas while the subscribed users have a single antenna. Extending the results to multiple 
receive antennas and using precoded spatial multiplexing is a subject of future investigation. We assume 
that the femtocells are distributed according to a homogeneous Poisson point process (PPP). We derive 
an expression for the probability of successful reception of limited feedback beamforming as a function 
of feedback delay, quantization size, the distribution and the density of the femtocell interferers in the 
network. The general mathematical framework is based on properties of the Chi-squared distributions 
and the homogeneous PPP. We use the probability of successful reception to compute a closed form 
expression for the maximum number of femtocells that can be deployed inside the macrocell as a 
function of the distance between the mobile user and the macro base station. 

We show that limited feedback beamforming outperforms no beamforming, in terms of average 
achievable rate. Thus we conclude that limited feedback beamforming is viable for transmission on 
the downlink of a femtocell network. To evaluate the achievable rate, we use the goodput metric p"3j, 



1 33 1, [34|. Goodput is defined as the number of bits successfully transmitted to the receiver per unit of 
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time. It takes into account the packet outage caused by the CSI mismatch between the transmitter and 
the receiver. When the probability of outage goes to zero, it is equal to the throughput. To combat the 
effect of packet outage, we implement a rate adaptation strategy based on rate backoff. To maximize the 
average goodput achieved, the transmitter, instead of transmitting at a rate corresponding to the estimated 
SIR, backs off the estimated SIR by a multiplicative factor, /3SIR, and modulates its signal based on the 
backoff SIR. The backoff factor (3 is chosen to maximize the average achievable goodput. Rate backoff is 



a suboptimal adaptive transmission strategy used in [35], [36|, to overcome channel estimation errors, in 
the context of centralized and distributed opportunistic scheduling for single input single output channels. 
In our analysis, we derive the achievable rate with rate backoff when the feedback delay is the main 
performance bottleneck, as well as when the mobile user is subject to cross tier interference. We show 
that rate backoff recovers the throughput scaling with SNR, and increases the average achievable rate, 
without requiring any coordination between the macrocell and the femtocells transmissions. 



C. Organization 

The rest of the paper is organized as follows. In Section [TTl we describe the limited feedback cellular 



network considered, and we review the assumptions used throughout the paper. Section III introduces 



the system goodput as a function of the cross tier interference and the feedback delay. In Section IV 



we derive a lower bound on the probability of successful reception at the mobile user as a function 
of the femtocell density, the distance from the base station and the feedback delay. Section [V] presents 
the rate backoff strategy suggested and derives the optimal backoff factor to mitigate the effect of 



interference and maximize the achievable goodput. Section VI presents numerical results that illustrate 
the performance of the limited feedback system and present the achievable sum rate using rate backoff 



at the transmitter. This is followed by concluding remarks in Section VII 



D. Notation 

Bold lowercase letters a are used to denote column vectors, bold uppercase letters A are used to 
denote matrices, non bold letters a are used to denote scalar values, and caligraphic letters A are used 
to denote sets or functions of sets. Using this notion, \a\ is the magnitude of a scalar, ||a|| is the vector 
2-norm, A* is the conjugate transpose, A T is the matrix transpose, [A]; m is the scalar entry of A in 
the £ th row and the k th column. We use E to denote expectation, a* to denote the metric a evaluated at 
the transmitter, and a to denote the average of a. 



II. System Model and Assumptions 

Consider a central macrocell Co overlayed with a tier of uncoordinated femtocells. The heterogeneous 
network considered is depicted in Figure [T] The macro base station B is located at the center of the 
macrocell, and is equipped with iV& antennas. Each femtocell is equipped with Nf antennas. The user 
terminals in the macrocell and the femtocells have a single receive antenna. We assume that each base 
station serves one active user at a time, using intra-cell orthogonal access. Both the macrocell and the 



femtocells employ a limited feedback beamforming system [37|, [38 1. The user terminal estimates the 
downlink CSI sequence using pilot symbols sent by its serving base station. We assume that the receiver 
estimates the channel perfectly. The CSI is quantized by means of a codebook and the quantization 
index is sent to the base station via a limited feedback channel. Since the macro base station and the 
femtocells are not synchronized and are assumed not to coordinate their transmissions, the quantization 
index at the user terminal only takes into account the downlink channel from its base station. The 
feedback channel is assumed to be error free, with a fixed feedback delay d. 

We use a stochastic geometry framework to model the distribution of the femtocells. The femtocells 
are assumed to be distributed according to a homogeneous spatial Poisson point process LTj with a 
fixed average density of A/ femtocells per meters squared. The average number of actively transmitting 
femtocells is hence Nof = A/C, where C = ixR 2 c is the area of the macrocell, R c is the radius of 
the macrocell. The PPP model, in addition to being analytically tractable, is justifiable by arguing that 
the femtocell transmissions are uncoordinated. It is the worst case among all stationary point processes 
due to its maximum entropy property. We assume that cross tier interference is the main source of 
interference. We treat the aggregate co-channel interference from neighboring macro base stations as 
noise, in the interest of evaluating the effect of the interference from the femtocells on the downlink of 
the macrocell. We focus on the transmission on the downlink of the macrocell to evaluate the feasibility 
of limited feedback beamforming in this scenario. We assume that the femtocell user is well shielded 
from interference on the downlink because of the indoor penetration loss, and the high received signal 
strength from its serving femtocell, as compared to the macrocell and the adjacent femtocells. 

The received signal at the cellular mobile user M is 

2/oN = s/^V^h* [n]f [n -d]s [n] + ^JJj ^ \^i\^sU n ] w i[ n } r i[ n ] +M n ], 0) 
where h [n] E C NbXl is the downlink channel between M and B , with zero mean, unit variance, i.i.d 
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complex Gaussian entries. The vector gj[n] G C NfXl is the downlink interference channel from the 
femtocell Fj to M . The vector f [n — d] G C^ 6 * 1 is the transmit beamforming vector at £? for M)- 
The vector Wj[n] G C NfXl is the transmit beamforming vector used by femtocell Fj to maximize the 
received power at its user terminal, for the transmitted signal The signals s and are such that 
E[|s | 2 ] < Pt, E[|rj| 2 ] < P t ., where P t and P ti are the maximum transmit power at the base station and 
the femtocells, respectively. p m and pf are the effect at the receiver of the transmit power, the carrier 
frequency and the wall partition loss, from the pathloss model. D is the distance between B and its 
mobile user M , and Dj is the distance between Fj and M . a rn > 2 and «/ > 2 denote, respectively, 
the pathloss exponent of the outdoor channel between B and M , and the indoor to outdoor channel 
between F and M . i; [n] ~ CAf(0,N ) is the additive white Gaussian noise at M with variance A^ . 

III. The Goodput Model with Limited Feedback 

In what follows, we derive the average goodput of limited feedback beamforming on the downlink of 
the macrocell network to study the viability of limited feedback in the heterogeneous network scenario. 
We then propose a rate backoff strategy to maximize the average achievable rate of the system. 

The macro cell user terminal M quantizes its channel h [n] by means of a unit norm vector codebook 
T = {f 1? f 2 , • • • , fjv}, where N = 2 B is the size of the codebook. The quantizer function Q chooses 
the beamforming vector i e that maximizes the SNR at M such that 

f e [n] = Q{h [n}} = argmax \h* [n]f k \ 2 , 1 < k < N. (2) 

f fe e.F 

The channel h [n] is then mapped to the index I n — £ which is fed back to B using B bits. 

The quantization process at M does not take into account the cross tier interference from the 
femtocells, as the macro base station and the femtocells are not synchronized, due to the unplanned 
deployment of the femtocells. M has no knowledge of the effective channels gjWj from F G II/. To 
modulate its signals, the base station estimates the signal to interference noise ratio (SINR) based on the 
delayed and quantized CSI as well as the statistical knowledge of the interference. This estimated SINR 
might be different from the actual SINR estimated at the receiver. When the transmitted rate exceeds the 
supported rate of the channel, a rate outage or a packet outage occurs. To account for the rate outage, 
we use the system goodput as the performance metric. The goodput gives the amount of information 
successfully received (without retransmission) at the user terminal. It provides a model to compute the 
effect of outage on the achievable rate, and permits the design of transmission techniques to combat 



the outage and consequently decrease the number of retransmissions needed. Goodput is defined as 

A[n] = R> - d]X (R*[n - d] < R[n]) , (3) 

where T(A) is the indicator function, which evaluates to 1 if the event A is true, and otherwise. 

While hybrid automatic repeat request (HARQ) has been recently used in cellular systems to combat 
the CSI mismatch problem, we assume in this paper that HARQ is not present. In systems employing 
HARQ, the signals received in error are retransmitted to maximize the probability of successful reception 
at the user terminal. This increases the packet delay and increases the receiver complexity. Moreover, 
the maximum number of HARQ retransmissions is fixed, and the long term goodput achieved with 



HARQ depends on the number of successfully received packets after the retransmissions Q39J, [ |40| . 

The rate supported by the channel R[n, d], for M , assuming Gaussian distributed transmit symbols 
so[n], is written as a function of the signal to interference ratio SIR[n, d], 



p m D-"^|hg[n]f [n-d]| 2 = |hg[n]fb["-d]| 2 
Pf 52 Fi en f IDin^lsrMwiMP Qd// 



R[n,d] =log 2 (l + SIR[n,d]) 

where SIR[n, d] 

We assume that the system is interference limited, and that the received signal power is much higher 
than the noise power. We thus ignore the thermal noise in the computations in this paper and mainly 
focus on the effect of cross tier interference on the achievable rate [|9j. Furthermore, the expressions 
using SIR[n, d] are more tractable analytically p2| . 

The interference function If = ^2 F . enf |Dj|~ Q/ |g*[n]w"j[n]| 2 is a shot noise process. The shot noise 
model is considered an efficient statistical model for predicting the fluctuations of the interference field, 
in the same way as the Rayleigh fading model is widely used to predict the fluctuations in the multipath 
channel, [41 1. The pathloss ratio Qd = — ^- am ls a function of the distance D, and the ratio of pf and 



p m . The transmitted rate is written in terms of SIR^n, d], the estimated SIR at B as, 

R t [n-d]=log 2 (l + SIR t [n,d]), (4) 

where SIR>,d] = p\h*[n - d]f [n - d]| 2 . 

We assume that the transmitter has an estimate of the average signal to interference power ratio 



The average goodput is expressed in terms of R*[n — d] and R[n, d] as 

A = E [R*[n - d]Z(R> - d] < R[n, d])] = E [R*[n - d]P [R^n - d] < R[n, d] | R*[n - d]]] 

/>oo 

= / R t (T)P[SIR > T | T]/ T (T)dT. (5) 
Jo 

The integration is taken over SIR^n, d], and /r(Y) denotes the probability density function (pdf) of 
SIR^n, d]. To evaluate A, we find a closed form expression for the probability of successful reception 
P [SIR > T], as a function of the delay, the femtocell density and the distance between M and B . 

IV. The Probability of Successful Reception 

The probability of successful reception is evaluated as the complementary cumulative distribution 
function (CCDF) of the desired channel power |hj!j[n]fo[Ti — d]| 2 given the interference shot noise // 



[SIR[n, d] > T] = 



|h5[n]f [n-d]| 5 



Qd \^i\- a, \snn}wi[r 
F,en f 



> T 



|hS[n]f [7i-d]| 2 > QdJ/T 



(6) 



where T denotes the SIR threshold. 

We use a Gauss Markov autoregressive model, shown in the literature to be reasonably accurate for 



small delays on the communications links [42|, [43|, to account for the temporal correlation in the 



system. The effective channel at the receiver, assuming a Gauss Markov model [44 1, is expressed as a 
weighted sum of two effective channels, 



h*[n]f [n - d] = (r}h* [n - d] + y/l - v 2 e*[n]) f [n - d], 



(7) 



where e[n] is a Gaussian vector with CjV(0, 1) entries, independent of h . The correlation coefficient r\ 
is determined using Clarkes isotropic scattering model as i] = j7o(27rd/dT s ), where f d is the maximum 
Doppler spread and T s is the symbol duration. The Doppler spread is fd = vf c /c, where v is the relative 
velocity of the transmitter-receiver pair, f c is the carrier frequency, and c is the speed of light. 

Lemma 1: The desired instantaneous channel power is written as a function of the first term of the 
Gauss Markov model h [n — d], for large values of i] corresponding to low mobility 



|h* [n]f [n-d]| 5 



^|hS[n-d]f [n-d]r. 



(8) 



Lemma 2: The desired effective channel power |h^[n]fo[n — d]| 2 can be expressed as the weighted 
sum of two independent random variables, Z = if ((1 — 5) X + Y), where X and Y are distributed as 
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X2(7Vf,-i)(l) and X2(l)» respectively. The CDF of the effective channel power is given by 

-i) + - xp (-A)gg_(A) , 

where Cl = (1 - 5) (l) N "~\ «i = 2rf (1 - 5), c 2 = (J)"'" 1 , and k 2 = 2rf. 

Proof: Using Gersho's conjecture on asymptotic quantization [ |45| , j |46j , it was shown in [ |47| that 
|h [n]| 2 sin 2 (Z(h [n], v)) has a Chi squared distribution X2{N b -i){$), and that |h [n]| 2 cos 2 (Z(h [n], v)) 
can be expressed as the sum of two independent Chi squared random variables, xi{N b -i){^-~ 8) and X2(l), 
respectively. Gersho's conjecture assumes that each quantization cell is a Voronoi region of a spherical 
cap with a surface area equal 2~ B of the total surface area of the unit sphere. This approximation has 
been shown to closely model the performance of codebook design techniques such as random vector 



quantization [47|, and is used in the literature to analyze the performance of limited feedback systems 



30J, [38 1, p7|. The distribution of the sum of two independent Chi squared random variables with 



2(Nb — 1) and 2 degrees of freedom, respectively, follows from [48 1. ■ 

Figure [2] compares the empirical CDF of the effective received power of |ho[n]fo[n — d]| 2 , with that 

of ?] 2 |hQ[n — d]fo[n — d]| 2 obtained in Lemma 2, for increasing velocities v 6 {10, 20, 30, 40, 50} km/h, 

for a 4 x 1 multiple input single output (MISO) downlink system, with a codebook size of B = 6. We 

observe that the approximation tightly follows the empirical CDF, with tolerable error up to moderate 

mobile speed of 50 km/h. 

Theorem 1: The probability of successful reception at a mobile cellular user, in the presence of cross 
tier interference, for a limited feedback MISO system, over temporally correlated channels, is given by 

P [SIR > T] = A 1 w 1 (T)e-^ l(T ) + A 2 e^ l(T) + c 2 e^ 2(T) , (10) 

where 

N b -2i-l w i-t-l , v /Tn x 2/ Q/ 

e E^-i)'-' n OH* - m = A ' c ' (^) • 

M = -c, g 4 <, „ 2 ,t, = x,c, (^2) 2 "" , »„d o, = (A) r _ A) . 

Proof: The proof is provided in Appendix |Aj ■ 
Using Theorem 1, we obtain an estimate of the maximum contention density A/ of femtocells in the 
network, given a probability of outage constraint e. 

Corollary 1: The maximum femtocell contention density A/(D) for which the probability of success- 
ful reception satisfies the maximum outage probability constraint P [SIR > T] > 1 — e, given T, is the 
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solution to 

A 1 w 1 (A / )e- Wl( ^' ) + A 2 e~ u ' lix f ) + c 2 e^ {x f ) > 1 - e. (12) 

In the high resolution regime, for small values of 5 corresponding to a large codebook size 2 B , A is 
given by 

Lambertw(- ^ ) + % + % 

A/ < ^ — A :, v , d3) 



-c 



where LambertW denotes the Lambert W function that solves the equation W / exp(W / ) = X for W as 
a function of X. 

Proof: The proof is provided in Appendix [B] ■ 
Estimating the maximum number of femtocells as a function of the distance of the mobile user from 
the macro base station allows the network designer to put a limit on the number of femtocells to be 
deployed in each area. The number of femtocells increases exponentially with the number of antennas 
at the macro base station and the quantization size for the limited feedback beamforming system. 

V. Rate Backoff 

When the estimated SIR 1 at B exceeds the estimated SIR at M , due to CSI mismatch, the transmitted 
rate exceeds the supported rate of the channel. This results in a packet outage. To decrease the probability 
of outage, and increase the probability of successful transmission, we implement rate backoff at the 
transmitter. Instead of transmitting at a rate corresponding to the estimated SIR*, the base station 
computes a backoff signal-to-interference ratio, SIR b (SIR 1 ), which backs off the estimated SIR 1 by 
a multiplicative factor (3 

SIR b = /3SIR\ (14) 

The backoff SIR* is chosen such that the average achievable goodput is maximized. Under the backoff 
model, the average goodput is written as 



A b =E 



log 2 (l + SIR b ) P[SIR b < SIR] 



(15) 



Setting the backoff factor, (3, to a value close to 1 provides a goodput rate close to that achieved without 
rate adaptation. A smaller (3 value, however, might be conservative and lead to a low average achievable 
rate. We are interested in computing the optimal value for /3 e [0 1] such that A b is maximized. The 
average goodput A b is a continuous and differentiable function in (3. Thus there exists an optimal backoff 
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factor (3* such that 

(3* = argmaxA b . (16) 



We first compute the optimal backoff factor due to feedback delay, without taking into consideration 
the femtocell interference. We then extend the analysis to include both delay and interference. 

A. Rate backoff with feedback delay 

In the absence of interference, the feedback delay is the only source of CSI mismatch at the transmitter. 
The rate supported by the channel is 

R N = log 2 (l + SIR)=log 2 (l + ' h °M f ^-^ 2 ) ! (17) 

l 



where, with no cross tier interference, the pathloss ratio is Q D = , 
The estimated SIR 1 at the transmitter is SIR* = l h c>- d ] f o["- d ll 2 _ 

In the absence of interference, the probability of successful reception at the receiver is given by the 
CCDF of the effective channel power 

-i)--p(-^)gg(^(^) • 

To maximize the average goodput with the suboptimal rate backoff function (3T, one needs to find a 
factor (3* such that the integral on the right hand side of 

/•OO 

A b =/ R*(/3T)P(SIR > /3T)/ T (T)dT (18) 
Jo 



is maximized. In ( 18), maximizing the argument of the integration results in maximizing the integration. 



We interchange the integration with the maximization, and the problem becomes finding (3* such that 

/?* = argmax[log 2 (l + /3T)P(SIR > /3T)] . (19) 





Theorem 2: The backoff factor (3* that maximizes the average goodput for a limited feedback beam- 
forming MISO system, with delay on the feedback channel, is the non-trivial solution of 

i _ s (Nb ' 1 K 



log 2 1 



(l + g) (kCACfl + 1 MP)) c 2 k 2 e^ 



gT 
Qd 



(20) 
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N b -2 i-l 



where f x {fi) = £ E 



(3T 



(i-£)\ V «i 
fci = — — , fc 2 = — — , and k-i 



JVi-2 i-l 



T 



i=l £=0 



/3T 



i-i 



Dlog(2)' 



For large codebook sizes, setting k 2 = «i, the backoff factor /3 can be computed as the largest root 
of the polynomial in (3 of order N b — 1 



dP?"- 1 ?) - C 2 P?»- 2 (/3) + C 3 P 3 Nb -\p) + C 4 /3 + C 5 = 



withp 1 ^- 1 (^)= ]T 2 



(21) 



i=0 «=0 



; ( — ] ,p^-\p) = P 3 ^- 1 (/3) = E 

£=0 



Ci = (1 - - C 2 = -(1 - 5), C 3 = T(l - 5)^- 2 , C 4 = (l + 2X) T-, and C 5 = 2<^X_ + X 

Proof: The proof is provided in Appendix [C] ■ 
The backoff factor /3* in the presence of feedback delay is a function of the distance D and the correlation 
coefficient r\. It also varies with the SIR threshold at the receiver and the number of feedback bits B, 
relative to the number of antennas at the base station, iV&. The backoff factor (3 admits an algebraic 
solution for Nf, < 6 by the Abel-Ruffini theorem. 



B. Rate backoff with cross tier interference and feedback delay 

In the presence of uncoordinated interference at the mobile user M , the probability of successful 



reception is given by (10). The backoff factor f3* computation follows similarly to the case of rate 



backoff with feedback delay. 

Theorem 3: The backoff factor (3* , that maximizes the average goodput for a limited feedback 
beamforming MISO system, in the presence of cross tier interference with maximum density A/(D) is 
the non-trivial solution of 



(Aiwi(/3T) + A 2 ) e-^W + C2e -^ T ) 



/3T 



log 2 



(l + [|e--^ T » ((* - A 2 )c^T) AiJUflT)) ~ ^»G8T)e 



■« a (/3T) 



(22) 



For large codebook sizes, the backoff factor (3* is computed as the 5/th square root of the roots of the 
polynomial in (3 5f of order 2, 



ArffultfT) + S f (c 2 - Ax(l + log(2)) + A 2 ) wi(/3T) 



A 2 +c 2 
log(2) : 



(23) 



such that (3* G [0 1]. 

Proof: The proof follows similarly to the proof of Theorem 2, with the expression of the probability 



of successful reception given by (10) 
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The backoff factor (3* is a decreasing function of the velocity of M and thus an increasing function 
of 77. This implies that as the velocity of the mobile user increases, the average achievable rate in the 
system, with rate backoff, decreases, and comes closer to the average achievable rate without backoff. 
It is further a decreasing function of the femtocell contention density A/, and the distance D between 
M and B . Thus, as the SIR at M increases, the average achievable rate with backoff increases. 

The backoff factor (3* is an increasing function of 5. Increasing the number of antennas, and fixing 
the feedback channel rate B, the average achievable rate with backoff increases relative to that without 
backoff. Similarly, increasing the rate on the feedback channel, for a fixed number of transmit antennas, 
results in a smaller backoff factor and a lower backoff rate. 



VI. Simulation Results and Discussion 

We consider a single macrocell of radius R c = 1km, overlayed with a tier of randomly distributed 
femtocells. The femtocells are distributed according to a homogeneous PPP with an average density of 
Xj femtocells per cell-site. We adopt a distance based pathloss model corresponding to the EVIT-2000 
channel model [49 ] for outdoor and indoor pathloss. For analytical simplicity, we do not consider random 
Log-normal shadow fading. We consider frequency flat Rayleigh fading with individual complex entries 
distributed as CJ\f(0, 1). We assume fixed wall partition losses corresponding to indoor-to-outdoor and 
outdoor-to-indoor propagation, equal to 5 dB. The outdoor and indoor-to-outdoor pathloss exponents 
are set to 3.8 and the carrier frequency is 2 GHz. The users are uniformly distributed inside each cell. 

We start by plotting the empirical CDF of the probability of outage at a mobile user with cross 
tier interference and delay. For a user average velocity of 20km/h and a target SIR threshold of 5 dB, 
Figure [3] compares the empirical CDF of the probability of outage versus SNR at the mobile user, to 
the expression of the probability of outage obtained using Theorem 1. The figure shows that the bound 
in Theorem 1 closely represents the probability of outage for all possible locations of the mobile user 
inside the macrocell. 



To get a practical sense of the number of femtocells obtained from Equation (12), we plot, in Figure 
4} the maximum number of transmitting femtocells Nf = TcR^Xf as a function of SNR at M , with a 



B 



10 percentile outage probability requirement, for increasing 5 = 2 ^s- 1 . The density A/ is averaged 
over 1000 uniformly distributed users in the macrocell, for an average user velocity of 20Km/h. As the 
number of feedback bits B increases, for example, for N b = Nf = A antennas, the number of femtocell 
interferers that can be allowed in the system without violating the probability of outage requirement 
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increases, this increase is more discernable for higher feedback rate, B > 6. Similarly, as the number 
of antennas at the base station N b increases, for a fixed number of feedback bits, A/ increases. 

Figures [5] and [6] illustrate the rate backoff concept applied at the transmitter to maximize the average 
goodput, respectively for delay and cross tier interference and delay. Figure [5] plots the ergodic goodput 
for a mobile user velocity of 20Km/h, N b = Nf = 2 and delay of 2 time frames. The goodput achieved 
when the probability of outage is equal to zero, i.e. the rate supported by the channel (the throughput) 
as a function of the SNR at the receiver is shown for comparison. We observe that the gap between the 
goodput and throughput is considerable. This gap increases with SNR as the rate of increase for the 
goodput depends on the probability of outage. Applying rate backoff at the transmitter recovers the rate 
of increase of the achievable rate with the SNR to that of the throughput. The gap between the backoff 
rate and that of the throughput remains because of the conservative transmitted rate, when the estimated 
SIR* falls short of the estimated SIR at M . We also plot the backoff goodput using the approximation 
in (22) for small § values. The approximation is shown to be accurate for N b = 2 antennas. It gives, as 
expected, a lower bound on the backoff factor (5* and the backoff goodput. 

Figure [6] shows the ergodic rate for Nof = 95. It compares the average goodput achieved with 
limited feedback beamforming, with and without rate backoff, to that achieved using open loop random 



beamforming [50|, for N b = Nf = 4 transmit antennas, and B = 5 feedback bits. The throughput is 
shown for comparison. We observe that limited feedback beamforming achieves an average rate gain 
of 5 dB over random beamforming for E[R[n]] = log 2 (1 + SINR). For the goodput A, this gain is on 
average 4 dB. This suggests that limited feedback beamforming, although prone to errors due to delay 
and quantization, is still a viable option for macrocell transmission use in heterogeneous networks MISO 
systems. Similarly to the system without interference, applying rate backoff at the transmitter recovers 
the rate of increase of the achievable rate to that of the throughput. The achievable rate with backoff is 
very close to that achieved with random beamforming, in the absence of outage. 

Figure [7] plots the backoff factor f3* as a function of the velocity of the mobile user and the number of 
feedback bits relative to the number of transmit antennas, 5. The figure shows that the backoff factor (3 is 
a decreasing function of the velocity v E [20 60] km/h. Rate backoff performs better for high temporal 
correlation coefficient -q, corresponding to pedestrian or low speed mobile users. This is reasonable 

B 

because the approximation in (8) assumes low mobility. j3* is an increasing function of 5 = 2 N t- i . 
Increasing the number of antennas for fixed feedback rate increases the backoff rate achieved. Similarly, 
decreasing the number of feedback bits, for a fixed number of antennas, decreases the backoff factor. 
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VII. Conclusion 

In this paper, we proposed a MISO limited feedback system for use on the downlink of a heterogeneous 
cellular system. We evaluated the effect of uncoordinated cross tier interference and feedback delay on 
the performance of the system. We derived a closed form expression for the probability of outage as 
a function of the distance from the mobile receiver to its base station, the velocity of the user and 
the femtocell density. We showed that, although prone to rate outage due to the inaccurate CSI at 
the transmitter, limited feedback MISO is a good candidate for use on the downlink of heterogeneous 
cellular networks. To maximize the average achievable rate, we proposed rate backoff techniques at 
the transmitter. We showed that rate backoff maximizes the goodput achieved, and hence increases the 
long term achievable rate in the system. Future work includes applying rate backoff to the second tier 
network considering both cross tier and same tier interference. 



Appendix A 
Proof of Theorem 1 
The probability of successful transmission is computed as 



>[SIR>T] = P[\h* [n]i [n-d}\ 2 > TQ D //, C ] 

sQdT 



c 2 exp - 



c 2 exp - 



ci exp I — 



«2 

£QdT 
sQdT 



d exp | — - — ) 2^ 2^ 



— ci exp 



N b -2 i-1 

EE 

i=\ 1=0 



sQ D T 



i=0 1=0 
N b -2 

E^ 



(i-iy. 



sQdT 



sQpT 
(i-£)\ V ki 

dP {If, a < s) 



i-i 



dP(//, c < s) 



(24a) 
(24b) 



where step (a) follows from conditioning on If c and substituting for the CCDF of |hQ[n]f [n — d]\ 2 . 

The integral in (24a) is expressed in terms of the Laplace transform of If c , £ Ifc (9) = E [e - ^ 6 *] = 
exp (-X f C f e 5 f) (from [22]), evaluated at 9 = ^ and 9 = respectively, for C f = 7v5 f T(5 f )T(l - 
6 f ), 

TQ D 



c 2 exp 



.sQpT 

K2 



. sQdT \ ^ 2 
d exp | — ) 2^ s 

i=0 



>(I ftC <s) = c 2 £ If:C 

V K 2 



/ N b -2 \ 

ci E <n c it. 

\ i=0 / 



c 2 exp -A/C/ 



V K 2 



c i 6 exp ~ f f 
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The integral in (24b) is evaluated as 



N b -2 i-1 



dP(i f , c < S )± Cl E 



i=l ^=0 



TQ D Y~* d 1 - 1 



d8 l 



N b -2 i-1 j- 



i=l f=0 



TQ D 



q5f — 



where (a) follows from the identity C [x k f(x)~\ = (—l) k F^ k \s) for the Laplace transform, with F^(s) 
representing the A;th derivative of F(s). Step (b) follows from evaluating the A;th derivative of the Laplace 



Transform C Ifc (6), [22 1 



d fc 



exp (-X f C f e 5 f) = - [XfC f n^} (S f - m)^- fe exp (-\ f C f 8 s f)] + &{\)C 2 f 9 25 s), 



and ignoring the second order terms. 

Combining the two parts of the integration, the probability of successful reception at M is 



'[SIR^T] = c 2 exp \-\fCf 



TQd 

k 2 



N b -2 i-1 



6< 



* E Et^v (~) [-A / C / n^o 1 (^ / -m)^-^exp(-A / C7 / ^)] 
i=l i=o ^ >' ^ 1 ' 

A 1 wi(T)e^ 1 ( T ' + A 2 e~ UJl{T) + c 2 e-^ T \ 



i-i 



(25) 



where A\, A 2 , wi(T), w 2 (T), and C/ are given in (11 ). 



Appendix B 
Proof of Corollary 1 

_ B 

For small values of 5 = 2 N b-*, corresponding to large codebook sizes, relative to the number of 
transmit antennas, we make the approximation Ki = k 2 = 2rf{\ — 5). This implies that wi(A/) = 
W 2(A/) = XfCf C^^*) S ■ The approximation gives a lower bound on the probability of success and 
consequently an upper bound on the probability of outage, c 2 e^ UJl( - XfS) < c 2 e~ u ' 2 ^ Xf \ The maximum 
femtocell contention density is derived from, 



^iw 1 (A / )e^ l(A ^ ) + A 2 e-^ l(x ^ + c 2 e^ l(A ^ > 1 - e 



C(A/)e c(A ^ > 



1 - e 



^2 , <=2 ' 



-A 1 e A i +A i 



where C(A/) = — o>i(A/) — ^ — It is given by the solution to the monotonically decreasing Lambert 



W function LambertW(A/) = C(\f)e c<yXf \ Consequently, A/ is given by (13) 
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Appendix C 
Proof of Theorem 2 

The expression for (3* follows from taking the derivative with respect to (3 of 

f3* = argmax[log 2 (l + /3T)P(SINR > /3T)] (26) 

i-e 



and setting the derivative equal to zero. Under the definitions = E<=i ^7 \ 77 

k\ = k 2 = an d k 3 = q d > tne optimization expression is written as given in Theorem 2. 
For large codebook sizes, we make the approximation k 2 = K\. This approximation allows us to 



remove the exponential terms from the expression in (20). Furthermore, we notice that 



i=0 £=0 v ; v 7 

' — ' V. Kl 

For Qd large, we have log 2 ( 1 + ^- ] ~ and — W w 1. The expression in (21 ) then follows by 

V / ^= D 1+ ^i7 



rearranging the terms in (20). 
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Fig. 2. The CDF of the effective received power |ho[n]fo[n — d]| 2 and of the approximation 77 2 |ho[n — d]frj[n — d]| 2 , for increasing 
velocities of the mobile user Mo v £ {10, 20, 30, 40, 50} km/h, for a 4 x 1 MISO downlink system with a codebook size of B = 6, and 
delay d = 2 time frames. 
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Fig. 3. The probability of outage P out = P[SIR < T] as a function of SNR or the distance D from M to B for N b = N f = 4, 
codebook size of B = 5, and mobile velocity of 20km/h. The analytical P out corresponds to that derived in Theorem 1. 




Fig. 4. The maximum number of femtocells interferers in a macrocell of area of C = ttIOOO 2 for a 4 x 1 MISO limited feedback system 
with increasing number of feedback bits B for Nb = Nf = 4, velocity v — 20 km/h, and delay d = 2 time frames. 
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Fig. 5. The average goodput for a MISO limited feedback system with AT/ = Nt, = 2, and B = 3 as a function of SNR with and 
without rate backoff. The femtocell density is assumed N a f = 95 femtocells per cell site. 




Fig. 6. The average rate as a function of SNR with and without rate backoff for a limited feedback system with TV/ = Nb = 4, and 
N f — 95 femtocells per cell site, and B = 5 bits on the feedback channel. 



24 




Fig. 7. The backoff factor (3* as a function of delta 5 £ [0.2 0.8], and velocity v £ [20 60] km/h, averaged over 100 users uniformly 
distributed inside the macrocell, with N f = 95 femtocells per cell site. 



